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Approximation  orders  of 
and 

approximation  maps  from 
local  principal  shift-invariant  spaces 

Amos  Ron 

1.  Introduction 

Let  S  be  a  function  space  consists  of  complex  (or  real)  valued  functions  defined  on  IR'*.  We 
say  that  S  is  shift-invariant  (SI.  for  short)  if  5  is  invariant  under  all  integer  translations  (referred 
to  hereafter  as  shifts),  i.o., 

(1.1 ;  Va  G  7Ld  (/  €  S  /(•  -  oi  G  S). 

In  this  paper  we  consider  SI  spaces  which  are  subspaces  of 

Lp  ;=  £p(l R'1). 

for  some  2  <  p  <  The  simplest  type  of  shift-invariant  spaces  is  the  PSI  space  ("P"  for 
■‘principar )  which  is  the  case  when  5  is  closed  (usually  in  the  underlying  p- norm,  hut  sometimes 
in  a  weaker  topology)  and  the  shifts  of  a  single  function  o  (=the  generator)  are  fundamental 
in  S.  Approximation  from  PSI  and  other  shift-invariant,  spares  is  pertinent  to  the  theory  ami 
applications  of  several  subareas  of  analysis,  and  in  particular  to  Multivariate  Splint'*.  Radial  Basis 
Approximation.  Wavelets  and  Sampling  Theory. 

In  many  actual  approximations,  the  SI  space  S  is  refined  to  yield  another  approximating 
space  Si,  with,  presumably,  better  approximation  properties.  The  standard  (known  as  stationary) 
refinement  is  by  scaling,  that  is.  S/,  is  obtained  by  dilating  the  functions  in  S: 

Sh  =  /r,.S  :=  {<Thf  :=  /(  /*)  :  /  €  S}. 

Sometimes  (cf.  [DR])  it  is  necessary  t)  refine  5  by  means  other  than  dilation. 

The  basic  way  for  measuring  the  approximation  “power"  of  5  is  via  the  tool  of  approximation 
orders.  Roughly  speaking,  the  collect),,,,  of  spaces  {S>,}/t>o  is  said  to  provide  approximation 
order  k  >  0.  if.  for  all  sufficiently  smooth  /, 

dist(/,  Sh)  —  0{hk). 

Here,  dist  is  measured  by  the  relevant  p-norm  or  one  of  its  relatives  (a  Sobolev  norm,  a  local  p-norm. 
etc.).  For  some  time,  the  analysis  of  approximation  orders  of  PSI  spaces  was  largely  dominated 
by  the  St  may -Fix  conditions.  [SF] .  These  conditions  assert  that,  if  S  is  generated  by  a  compactly 
supported  o.  if  df  0)  -£  0.  and  if  the  scale  {Sh  K  is  stationary,  then  the  approximation  orders  provided 
by  { St, }/,  are  determined  by  the  order  of  the  zero  S  has  at  each  of  27rZZri\0.  The  standard  method  for 
converting  the  information  about  these  zeros  into  approximation  order  results  is  the  polynomial 
reproduction  /  quasi-interpolation  argument;  (cf.  the  book  [C],  the  survey  [Bj.  and  the 
references  therein).  However,  several  important  PSI  spaces  that  were  introduced  and  Studied  in 
recent  years  do  not  satisfy  the  requirements  imposed  above  on  the  PSI  space.  One  difficulty  arises 
in  the  area  of  radial  basis  functions,  since  there  the  typical  generator  <j>  is  not  compactly  supported. 
A  totally  different  difficulty  arises  in  the  area  of  box  splines:  while  the  box  spline  is  compactly 
supported,  its  corresponding  {S/,}/,  is  not  a  stationary  one  (unless  the  generator  <3  is  a  polynomial 
box  spline).  The  attempts  to  cover  those  cases  by  generalized  quasi-interpolation  arguments  led 
to  some  remarkable  achievements,  but  did  not.  solve  the  problem  in  its  entirety.  In  retrospect,  it 
seems  that  the  quasi-interpolation  approach  fails  to  realize  the  approximation  order  of  general  PSI 
spaces,  a  fortiori  of  general  SI  spaces. 
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New  approaches  that  circumvent  quasi-interpolation  were  recently  developed  in  (BR2j  i p  =  x 
and  [BDRlj  [p  =  2).  While  the  two  methods  differ  in  the  approximation  scheme  they  emplov  a: 
well  as  in  their  error  analysis,  they  both  perform  the  entire  error  analysis  on  the  Fourier  domain 
hence,  in  turn,  avoid  the  imposition  of  decay  rates  on  the  generator  o.  Also,  importantlv.  pot! 
eharae'erizations  do  not  require  {S.ij}/;  to  he  stationary,  and.,  furthermore.  [BDRlj  even  dispeusi 
with  the  regularity  assumption  oil))  /  0.  Extensive  discussions  of  the  various  applications  of  tin 
results  of  [BR2]  and  [BDR1]  to  radial  function  approximation  can  he  found  in  ,3  of  [BR 2]  am 
[R3j  respectively.  Further  discu»sion  of  the  literature  can  be  found  in  §3.2. 

In  the  present  paper,  we  revisit  PSI  spaces  generated  by  compactly  supported  functions  isuei 
spaces  wore  coined  local  in  [BDR2]).  Primarily,  we  aim  at  establishing  results  on  the  approximate! 
order  of  bi'x  sphnr  spaces,  ami  providing  explicit  approximation  schemes  that  realize  that  order 
The  main  findings  here  with  respect  to  box  spline  spaces  are  as  follows: 

a : :  Wo  determine  exactly  (in  Theorem  3.-1)  the  approximation  order  of  box  spline  spaces  in 
the  Ly norm.  In  no  norm  has  such  result  been  known  before,  although.,  for  p  =  x.  [BR2]  came 
very  close  to  the  mark.  We  also  show  that  the  saturation  class  is  trivial  ief.  Theorem  3.7  s 
■;bi:  We  present  sin  §3.3)  an  explicit  approximation  scheme  in  the  form 


z 


a)[.ff)  (a). 


tor  a  suitablv  da.st  ;;  convolution  operator  f.  We  th  ♦  scheme  realizes  th.-  £>-approxh....;;«, 
order,  and.  as  a  matter  of  fart,  maintain  the  same  approximation  order  in  Lv- norms.  2  <  p  <  x 
.  Theorem  3.13  .  The  results  also  apply  to  the  I  x -ease.  if  a  mild  smoothness  condition  is  imposed 
on  o  -  O  €  A.  with  .4  the  Wiener  algebra). 

•  e i:  We  (ievelop  a  new  error  analysis  method  which  differs  from  the  two  of  [BR2]  and  [BDRlj. 
With  the  aid  of  that  approach,  we  show  that  for  2  <  p  <  >  twit  It  some  additional  smoothness 
conditions,  requirt'd  of  the  box  spline  o.  in  the  />  =  x  case)  the  approximation  scheme  used  provides 
simultaneous  approximation  to  functions  and  their  derivatives  (Theorem  3.17). 

Most  of  the  analysis  that  is  developed  to  deal  with  box  spline  spaces  was  found  to  apply 
a!-o  to  general  local  PSI  spaces.  Due  to  that  reason,  we  first  present,  in  §2.  results  that  concern 
approximation  orders  of  local  PSI  spaces,  and  only  then  apply  'hese  results,  in  Jj3.  to  the  box  spline 
case.  Most  of  the  proofs  are  collected  in  §-l. 

Soon  after  rumpiering  the  present  paper.  I  received  a  preprint,  of  Kyriazis"  paper  [K].  The 
paper  [Kj  deals  with  similar  problems  to  those  considered  in  the  present  paper,  and  to  a  limited 
ex*en?  employs  similiar  analysis  :  primarily,  both  papers  invoke  the  approximation  scheme  of  [Bl<2j. 
ami  also  both  paper  apply  the  Hausdortf- Young  theorem).  However,  the  foci  of  the  two  papers 
seem  p>  be  quite  different,  hence  they  seem  to  be  complement  ary  one  to  the  other.  [KJ  focuses  on 
stationary  refinements,  and  attempts  to  impose  as  mild  as  possible  decay  condition  of  the  generator 
O.  Therefore,  a  substantial  effort  is  devoted  there  to  proving  that  various  constructed  approximates 
are  indeed  taken  from  the  underlying  PSI  space.  The  paper  [K]  suceeds  in  providing  satisfactory 
results  also  to  £  norms  where  p  <  2.  its  class  of  smooth  functions  is  less  restrictive  than  here,  and 
if  established  Strang-Fix-like  conditions  similar  to  those  derived  in  [BDRlj.  The  present  paper 
focu-es  on  box  spline  space,  hence  avoids  on  the  one  hand  questions  of  decay  rates  of  the  generator, 
hut  needs  on  the  other  hand  to  deal  with  non-stationary  refinements.  Also,  we  present  here  results 
mi  simultaneous  approximation  to  functions  and  derivatives,  a  topic  which  was  not  considered  in 
[Kj.  and.  finally,  our  results  apply  to  p  —  \.  as  well  ( [K]  relies  on  the  Miehiin  multiplier  theorem, 
hence  deals  with  p  <  oc). 

Notations:  we  have  already  used  the  notation  .jy,  for  the  dilation  operator 

:/*—/( -/*) • 


Norms  of  vectors  x  6  IR'*  are  denoted  by  |xjp,  namely. 

d 

kip  :=(S>jn1/P’ 

3  =  1 

with  the  default  notation  |x|  :=  |x|_>.  The  function 

x  ~  |xj, 

which  is  used  extensively  in  the  paper,  is  denoted  by  ^ t ht>  essentially  self-understood  notation) 

H 

For  8  6  €‘l,  the  notation  eg  stands  for  the  exponential  function 

eg  :  at  - 

Unless  otherwise  stated,  all  domains  of  functions  in  this  paper  are  taken  to  be  IR'\  Thus.  LP  ~ 
LP  UV‘i.  S'  =  S’(JRJ)  <the  space  of  all  d-dimensional  complex-valued  tempered  distributions), 
liy  =  U’p  (1R'<)  (  the  Sobolev  space  of  all  functions  whose  derivatives  up  to  order  k  art'  in  Lp).  etc. 
We  also  abbreviate 

li/iir  :=  ii/lkp- 

2.  Approximation  from  local  PSI  spaces 

Our  model  is  as  follows.  W'c  are  given  an  indexed  set,  <f>  :=  {o,,  }A  C  Lp.  The  locality  assumption 
usually  means  that  each  0/t  is  supported  in  some  bounded,  ft-independent  domain  S2.  but,  while 
such  art  assumption  holds  indeed  in  the  box  spline  case,  we  do  not,  need  it  here.  We  only  assume 
that  each  6 h  is  compactly  supported.  Regardless  of  the  value  of  p.  we  define,  for  any  compactly 
stipported  o.  the  PSI  space  Sio)  to  be  the  infinite  span  of  the  shifts  of  <p: 

S{0)  :=  {  ®('  “ 

aeZ'.d 

The  convergence  of  the  infinite  sums  can  be  taken  pointwise.  since  the  sum  is  actually  finite  on 
compact  domains.  Xo  a-priori  growth  condition  is  imposed  on  the  coefficients  f ''(<'*)  h>e2U-  Al¬ 
though  this  definition  slightly  deviates  from  the  one  given  in  the  introduction  (our  space  is  not 
a  snbspace  of  Lp.  nor  S(o)  O  Lp  need  be  closed),  that  difference  would  not  matter  in  subsequent 
discussions. 

The  scale  of  spaces  {5h}h  is  obtained  by  dilating  the  PSI  spaces  S(Sh)' 

Sh  :=K/  =  /(/ft):  /€S(< M), 

and  the  approximation  orders  provided  by  {©/,},,  is  concerned  with  the  rate  of  decay  of 
(2.1 )  distp(/,  Sh)  :=  inf{||/  -  si!,,  :  a  €  5*} 

as  ft  —  0.  More  precisely,  we  say  that  {0/;}h  provides  approximation  order  k,  if  for  every  /  in  some 
smoothness  space  Vp  ^  and  small  enough  ft, 

dist.p(/,  Sh)  <  const  ftfc||/||p.fcT 


with  !! S'\p.k  some  norm  of  /.  The  scale  {S*}*  is  stationary  if  oh  —  o.v  for  all  h,h’ .  in  such  a 
case.  {Sh\h  are  all  dilates  of  one  basic  PSI  space  S{o). 

The  space  I'pit  of  "smooth  enough"  test  functions  is  defined  as  follows.  For  two  conjugate 
exponents  1  <(/</)<  oc.  and  k  >  0. 

12.21  Ip.*  :=  {/  :  ||/!U  :=  lid  +  \-})kf\\W^)  <  *)■ 

Note  that,  for  an  integer  k.  the  Hausdorff- Young  Theorem  implies  that  \pk  is  continuously  embed¬ 
ded  into  the  Sobolev  space  U  *. 

Given  f  €  we  seek  an  approximate  for  /  from  S*.  Since  S/,  is  the  /i-dilate  of  S(or,)-  we 
can  define  the  approximan*  for  /  in  terms  of  an  element  A-nif  >  €  S(o>, i.e..  approximate  /  In 
fii'-’W./)'.  i/'  is.  necessarily,  a  'possibly  infinite)  linear  combination  of  the  shifts  of  o,>,.  We 

obtain  the  coefficients  in  this  combination  as  the  restriction  to  7Z'1  of  a  continuous  1  in  fact,  entire; 
function  Jkif  '■  In  summary,  we  approximate  /  by  a >,.4 »,(/),  where 

i2.:P  Ah{f)  =  ]T  <)*(•- a)  J*(/)(a). 

Tims,  tin*  particular  details  of  our  approximation  scheme  rely  on  the  choice  of  the  maps 

/  -  jh[f). 

As  it  turns  out.  the  results  below  on  approximation  orders  requite  four  conditions  of  the  maps 
y,  ”.td  any  collection  that  satisfies  these  four  properties  will  do  here.  Three  of  these  conditions 
are  independent  of  the  specific  approximation  order  we  are  after,  and  are  listed  now. 

(2.4)  Conditions  required  from  the  maps  { ■/*. } ,h : 

(a)  Each  Jh  is  a  dilation  followed  by  convolution,  that  is 

<2.T)  •/*(/)  -  (7),^/)v. 

Ail  Each  Tit  is  a  function  supported  in  some  /i-independeut  origin-neighborhood  B  C  [ —  tt  . . 
id  For  some  hi j  >  0.  1„  are  uniformly  bounded  on  B  (lienee  on  IR'1). 

Note  that  we  do  not  impose  smoothness  conditions  on  71, .  and  therefore  ./i,  need  not  to  map 
Lp  into  itself.  However.  {a>,Jh}h  are  uniformly  bounded  etidomorpiiisms  on  each  Vp^. 

Theorem  2.6.  Let  2  <  p  <  oc  be  given  and  let  q  be  its  conjugate  exponent.  For  k  >  0.  let  \ k  be 
as  in  (2.2).  Let  {oh}h  be  a  family  of  compactly  supported  functions,  and  B  an  origin-neighborhood. 
Assume  that  the  collection  of  sequences 

rt>k,h  ■  (2~2Z,i\0)  3  3  c—  ]|(/i  4-  |  |)  +  *3)ll/JOC(/j).  h  <  ho. 

lies  in  Lq(2tt2Z'l\0)  and  is  hounded  there.  Suppose  that  {.Jh  satisfy  conditions  fa),  (b)  and  ( c ) 

of  (2.4)  (with  respect,  to  the  present.  B).  and.  in  addition, 

Id;  sup  || (h  +  J-|)“fc(l  -  <l>aTh)\\[,vtB)  <  oc. 

h</i0 

Let  Afl  be  defined  by  (2.3).  Then: 

11/  ~  0h(Ah(f)) <  const  hk\\f\\vk.  /  6  VP,k- 


Proof.  of.  §4.1.  4 

We  remark  that  the  proof  of  Theorem  2.6  provides  the  following  bound  on  the  error: 

11/  -  °>i-'U(/)!!;.i)(!RJ)  < 


t--")  const  (l!T/lj]i,K(S)!|m/c./,!|/i;(2!rz‘‘\o)  +  IK*  +  I'D  ‘(1  -  +  «(1/). 

with  the  o(ll  expression  always  bounded  by  1.  decays  to  0  with  h.  and  otherwise  depends  onlv  on 
/,  p.  k  and  B.  and  with  const  depending  only  on  p.  Therefore,  assuming  that  B  is  fixed,  one  might 
try  to  choose  -4  such  that  the  sum  +  i!i,!  +  1  -  Oj, 7V ) M /. .  'ft  is 

minimized. 

A  natural  choice  for  JhL  is  given  by  7),  :=  \joi ,,  with  \  a  characteristic  function  of  some 
0- neighborhood  B.  In  this  case  {h  4-  j-|)— *(  1  -  CO,  74, )  vanishes  on  B.  hence  condition  fdi  in  the 
theorem  trivially  holds.  The  only  condition  that  needs  to  be  checked  then  is  condition  (<:).  viz.,  the 
uniform  boundedness  of  which  amounts  to  the  uniform  boundedness  away  of  zero  of  l/Oh- 

Hence  we  have: 

Corollary  2.8.  Assume  that  the  Fourier  transforms  of  the  family  'h  :=  {<:>!,_} h<k,-,  of  compactly 
supported  functions  are  uniformly  bounded  away  from  0  on  sonic  origin-neighborhood  B.  Then, 
for  every  2  <  p  <  x.  every  k  >  0  and  h  <  hv. 

distpC f.Sh)  <  coP"t/ifc||/!|p.fc(||m*.h|!,i|(2.TZ<\o)  +o(I)).  V/  €  \’pk. 

with  {rnic.h}  defined  as  in  Theorem  2.6.  with  const  independent  of  k  and  f.  the  o(l,i  expression 
bounded  by  I.  and  with  q  the  conjugate  of  p.  Hence.  {ok}h  provides  approximation  order  no 
smaller  than  k  whenever  }/,  are  uniformly  bounded  for  sufficiently  small  h. 

In  order  for  this  corollary  to  be  useful  in  the  derivation  of  approximation  orders  we  need  to 
find  conditions  which  guarantee  the  boundedness  of  the  sequences  {m/c.a}/,-  In  the  next  section,  we 
will  see  how  this  is  done  in  the  case  of  box  splines.  At  present,  we  note  that  the  essential  part  in  the 
boundedness  assumption  of  {nik,h}h  is  the  pointwise  boundedness,  that.  is.  for  every  3  €  2~7Zd\0 
the  function 

A  -  \\(h  +  \-\rkOh(-±3)\\L^B),  h<  h0 

should  be  bounded,  and  the  bound  should  be  uniform  in  3.  The  fact  that  we  assume  more  than 
that  pointwise  boundedness  in  the  theorem,  is  due  to  the  technical  details  of  the  proof,  and.  in 
most  practical  examples  is  translated  to  smoothness  conditions  on  {Oh}s-  Note  also  that  in  the 
stationary  case.  oh  does  not  change  with  ft,  and  the  pointwise  boundedness  condition  thus  can  be 
easily  seen  to  be  equivalent  to  d  having  a  fc-order  zero  at  each  of  3  €  2T2Zd\0. 

Under  additional  smoothness  conditions  on  {c4}a,  the  approximants  {At,(f)}>,  can  be  shown 
to  approximate  /  in  Sobolev  norms  as  well.  We  mention  that  such  results  (concerning  simultaneous 
approximation  from  SI  spaces)  are  a  rarity,  especially  since  there  is  no  standard  way  to  derive  them 
from  quasi-interpolation  arguments.  The  most  notable  exception  is  [SF]  that  states  such  results  in 
the  L ■)-  and  norm  (for  the  stationary  case),  and  proves  the  /^-statement. 

Theorem  2.9.  Adopting  the  notations  and  assumptions  of  Theorem  2.6,  assume,  in  addition,  that 
for  sortie  positive  integer  r  <  k,  the  sequences 

mrkM  :  (2aZZrf\0)  33- mil h  +  |-|)-‘^(-  +•  0)\\lx{B),  h  <  A0 
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lie  in  tg[2~7Z.,i\ 0)  and  are  uniformly  bounded  there.  Then. 


||/  -  crh(Ah(f))\\w;  <  constp  r!!/;|pjtA* 


Proof.  of.  §4.1. 


3.  Approximation  from  box  spline  spaces 

3.1.  Statrnent  of  the  problem  and  its  /..--solution 

To  define  a  box  spline.  we  let  E  be  a  rational  matrix  of  d  rows  which  is  also  considered  as  *  h« 
multiset  of  its  columns  with  each  column  /  €  E  (referred  to  sometimes  as  "a  direction";  b 

assumed  to  be  a  non-zero  vector.  The  matrix  E  is  augmented  by  a  row  vector  A  =  A=  €  C~.  am 
the  resulted  matrix,  denoted  by  (E.A).  is  used  to  define  the  box  spline  .1/  :  =  .l/r  \  whose  Fouriet 
transform  is 

_  _  -it  *  _  1 

1 3. 1  M (~)  =  17  - - - .  w  6  1R,(. 

_  Af  -  if  ■  ^ 
ten  -  s 

In  general  .1/  is  a  compactly  supported  measure  defined  on  IR'f .  but  upon  assuming  that 
(3.2;  rank  E  =  </ 

las  we  always  do  hereafter),  the  box  spline  is  a  bounded,  compactly  supported  pierewise-exponentiui- 
poivnomiai  function  supported  itt  the  zonotope 

(3.3)  :=  iY  t(Z  :  U  €  [0  . .  lj }. 


Polynomial  box  uplines  correspond  to  the  choice  A  =  0.  Exponential  B-spl ines  are  obtained  when 
d=l  and  /  =  1.  all  (;.  Tensor  splines  are  obtained  whenever  all  the  directions  are  standard  unit 
vectors.  The  box  spline  is  positive  in  the  interior  of  Zz  whenever  A  is  real-valued. 

We  now  define  the  scale  {Sy}y  of  box  spline  spaces.  For  this,  we  fix  M  (i.e..  fix  E  and  A),  and 
define,  for  a  given  refinement  parameter  h  >  0,  the  box  spline  My  as 

Mh  ■=  Mz,h\. 

The  rest  of  the  definition  is  as  in  the  introduction,  i.e..  Sy  :=  <7yS(My).  Our  space  of  "test 
functions"  remains  the  space  \'p  y  defined  in  (2.2). 

Since  the  ladder  of  spaces  is  determined  as  soon  as  M  is  chosen  (in  affect,  as  soon  as 

E  and  A  are  chosen),  we  refer  to  the  relevant  approximation  orders  as  provided  by  M.  rather  than 
"provided  by  {A//}".  Note  that  each  Sy  is  “spanned"  by  the  hZZ  -shifts  of  the  dilated  function 
oy.Uy.  Furthermore,  in  case  A  =  0.  My  =  A/,  all  h.  and  the  scale  {.S\}A  becomes  stationary. 
The  motivation  behind  the  particular  definition  of  Sy  in  the  non-stationary  case,  is  that,  while  Sy 
becomes  invariant,  under  finer  and  finer  shifts  as  h  — -  0,  the  functions  in  Sy  are  always  piecewise 
in  some  finite-dimensional  E-independent  space  H  (of  exponential-polynomials). 

We  have  seen  in  the  last  section  that  approximation  orders  from  PSI  spaces  can  be  understood 
in  terms  of  the  behaviour  of  the  various  generators  around  2-zZ'i\0-  In  t.he  box  spline  case,  however, 
•inch  results  cannot  be  considered  as  satisfactory:  the  itntnediatelv  available  information  on  the  box 
spline  is  the  matrix  (E.A).  and  therefore  we  wish  to  characterize  the  approximation  order  of  box 
spline  spaces  in  those  terms:  that  is.  given  k  >  0.  we  need  to  find  all  (E,  A)  whose  corresponding 
box  spline  \fz.x  provides  an  approximation  order  k. 

For  the  L^-norm,  we  provide  in  this  paper  the  following  complete  answer  to  the  above  problem. 


Theorem  3.4.  The  L2-approximation  order  provided  by  the  box  spline  Ms, a  is  the  number 


(3.5)  k'(3)  :=  min{#Kp  :  (3  €  27r2Zd\0}, 
with 

(3.6)  Kp  :=  Kp{ 3)  :=  {£  G  E  :  £  •  0  G  2ttZ\0}. 

In  particular,  the  approximation  order  is  independent  of  A. 

Proof.  cf.  §4.3.  4 

Note  that  for  p  =  2  and  integer  k,  V2,k  =  IV2,  and  hence  the  above  stated  approximation 
orders  apply  to  the  entire  Sobolev  space. 

We  will  also  show  that  the  saturation  class  associated  with  the  above  problem  is  trivial.  Pre¬ 
cisely,  we  have: 

Theorem  3.7.  Let  {Sh}h  be  the  box  spline  space  scale  associated  with  a  box  spline  Ms, a-  Let 
k  :=  k'(3),  and  let  /  G  IF* \{0} .  Then,  for  every  sequence  hj  — »  0, 

j— °° 


dist 2(f,Sh])  #  o(h*). 


Proof.  cf.  §4.3.  4 

The  definition  of  k'(3)  is  entirely  in  terms  of  the  matrix  3  (i.e.,  does  not  require  any  information 
on  the  underlying  box  spline  A/s,  a),  and,  moreover  k'(3)  can  be  computed  by  a  finite  algorithm. 

In  the  important  special  case  when  3  is  an  integer  matrix  Theorem  3.4  implies  the  following 
(known)  result: 

Corollary  3.8.  Assume  that  3  is  an  Integer  matrix.  Then  the  L2-approximation  order  provided 
by  the  box  spline  A/e,a  is 

(3.9)  Jfc(E)  :=  mint#*  :  *  C  H,  rank(E\*)  <  d }. 

Here  and  hereafter,  XcE  means  that  X  is  obtained  from  E  by  the  deletion  of  some  columns, 
and  #*  is  the  number  of  columns  in  *. 

Proof.  In  view  of  Theorem  3.4,  it  suffices  to  show  that,  for  an  integer  matrix  E,  k'(3)  — 
k(3).  Let  £  G  E.  Since  £  is  integer,  £  •  0  G  27rZZ  for  every  0  G  27rZ5d\0,  hence  Kp  of  (3.6)  can 
equivalently  be  defined  here  as 

(3.10)  ^  =  {(eS:(-^0}. 

Thus,  ( 3\Kp)TP  —  0  and  hence  rank(E\A^)  <  d.  This  show  that  k(3)  <  #Kp,  and  consequently 
Jfc(E)  <  k'(3).  The  reverse  inequality  does  not  require  the  integrality  of  E,  but  only  its  rationality: 
assume  that,  for  some  X  C  E,  rank(E\*)  <  d.  Since  E  is  a  rational  matrix,  (E\*)  is  rank-deficient 
if  and  only  if  there  exists  a  non-zero  integer  vector  a  perpendicular  to  all  £  G  (E\* ).  In  view  of 
(3.6),  we  have  K2na  C  X,  and  hence  A:'(E)  <  #*2*0  <  #*•  It  follows  that  k'{ 3)  <  A;(E).  4k 


7 


3.2.  Further  literature  discussion 

Now  that  the  two  numbers  fc(E)  and  k'( E)  are  introduced  and  their  connection  to  approxima¬ 
tion  orders  is  revealed,  we  are  able  to  discuss  the  history  of  the  problem  in  further  detail.  In  this 
regard,  it  seems  instructive  to  separate  the  discussion  of  the  polynomial  box  spline  case  (A  =  0)  from 
the  general  exponential  case.  As  mentioned  before,  the  problem  of  the  former  case  is  stationary,  i.c., 
the  spaces  {S/J  are  all  obtained  from  the  original  space  S(M)  by  dilation,  and  there  is  a  variety  of 
papers  (including,  but  not  restricted  to,  [SF],  [DM1,2],  [BJ],  and  [JL])  which  treat  such  setting  for 
a  general  compactly  supported  0,  and  links  restrictive  (hence  stronger)  notions  of  approximation 
order  (known  as  “controlled”  and  “local”)  to  the  polynomials  in  S(M).  Further,  more  recently, 
it  was  shown  in  [R2]  (Lx- norm)  and  [BDR1]  (I2- norm)  that  whenever  0  is  compactly  supported 
and  0(0)  7^  0  (which  is  certainly  the  case  for  a  polynomial  box  spline  0  =  M)  the  polynomials 
in  S(M)  characterize  the  unconstrained  approximation  order  (i.e.,  the  one  defined  and  analyzed 
in  the  present  paper).  Thus,  at  least  in  essence,  the  characterization  of  the  approximation  order 
provided  by  a  polynomial  box  spline  amounts  to  the  identification  of  the  polynomials  in  S(M). 
These  polynomials  were  characterized  by  de  Boor  and  DeVore  in  [BD]  for  the  three-directional 
polynomial  box  spline  (d  =  2,  A  =  0,  £r  6  {(1, 0),  (0, 1),  (1, 1)},  €  E)  which  were  also  introduced 

there.  Polynomial  box  splines  associated  with  a  general  integer  E  were  introduced  and  studied 
by  fie  Boor  and  Hollig  in  [BH] ,  with  the  identification  of  the  underlying  polynomial  space  being 
among  the  highlights  of  that  paper.  The  abstract  argument  provided  in  [BH]  for  the  conversion 
of  the  knowledge  on  the  polynomials  into  lower  bounds  on  the  approximation  order  has  become 
a  standard  tool  since  then.  Another  proof  of  that  result  is  included  in  the  subsequent  work  of 
Dahinen  and  Micchelli,  [DM1].  The  characterization  of  the  approximation  order  of  a  polynomial 
box  spline  associated  with  a  general  E  was  only  recently  established  in  [RS],  where,  again,  the  main 
result  is  concerned  with  the  identification  of  the  polynomials  in  S{M). 

Exponential  box  splines  were  introduced  in  [Rl],  and  that  paper  also  contained  the  first  result 
on  their  approximation  order  (showing  that  for  A  C  1R  and  an  integer  E,  the  approximation  order 
in  the  Loo-norm  is  at  least  1).  The  first  comprehensive  discussion  of  approximation  orders  for 
exponential  box  splines  is  found  in  [DR],  where,  for  general  A  but  integer  E,  fc(E)  was  proved  to  be 
a  lower  bound  on  the  L^-approximation  order.  (It  was  further  shown  there  that  the  exponential 
reproduction  argument  cannot  provide  better  bounds).  The  extension  of  these  results  to  p  <  oo 
was  done  by  Lei  and  Jia  in  [LJ],  where,  in  addition,  the  local  structure  of  the  spline  space  was  used 
to  provide  matching  upper  bounds  (Thus,  Corollary  3.8  is  [LJ]’s). 

As  already  alluded  to  before,  all  the  aforementioned  results  employed  the  quasi-interpolation 
argument:  first,  the  space  H (A/)  of  all  exponential-polynomials  (polynomials,  if  A  =  0)  in  H/,5/,  is 
computed  (either  explicitly  or  as  the  kernel  of  explicit  differential  operators)  and  then  the  approx¬ 
imation  power  of  H(A/)  around  the  origin  is  studied.  This  local  approximation  order  of  H(A/)  is 
them  converted  to  lower  bounds  on  the  approximation  orders  via  the  quasi-interpolation  argument 
(for  this  argument,  in  the  exponential  case,  see  [DR],  [LJ],  [R2],  [BR1]  and  [CW]).  However,  in 
contrast  with  the  stationary  case,  there  is  no  general  theory  that  can  be  applied  to  the  exponential 
box  spline  to  show  that  the  lower  bounds  obtained  by  quasi-interpolation  are  the  best  approxima¬ 
tion  orders.  Indeed,  we  draw  (in  §3.3)  an  example  of  a  box  spline  M  (necessarily  with  a  non-integer 
E  and  a  non-zero  A)  such  that  k'(S)  =  1  while  the  corresponding  H(A/)  is  trivial.  We  stress, 
however,  that  examples  of  this  type  are  the  exception  rather  than  the  rule. 

The  only  reference  that  we  are  aware  of  which  treats  the  approximation  order  provided  by  M-=,x 
for  a  rational  E  and  general  A  is  [BR2],  where  approximation  in  oo-norm  is  considered.  The  ap¬ 
proach  in  that  paper  is  based  on  the  theory  for  Loo-approximation  orders  developed  there,  a  theory 
which  indeed  circumvents  quasi-interpolation.  Theorem  3.13  is  similar  to  its  Loo-counterpart  from 
[BR2],  with  one  important  difference:  while  the  results  of  [BR2]  require  some  minimal  smoothness 
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conditions  of  the  basis  function  under  consideration  (here  the  box  spline  M),  hence  exclude  box 
splines  of  low  smoothness,  no  such  exclusion  exists  in  I2-analogous  results  from  [BDR1].  We  will 
elaborate  on  this  point  in  the  next  subsection,  since  Theorem  3.13  and  its  proof  provide  a  better 
understanding  of  the  nature  of  the  smoothness  restriction  on  M  which  was  required  in  [BR2]. 


3.3.  Approximation  orders  of  box  spline  spaces  in  Lp,  p  >  2 

In  order  to  derive  approximation  orders  for  box  spline  spaces  we  invoke  Theorems  2.6  and  2.9. 
For  that  task  we  need  first  to -find,  for  the  given  box  spline  scale  {<f>h  =  A/*}*,  corresponding  maps 
(.//,}/,  that  satisfy  the  four  requirements  specified  in  Theorem  2.6.  Upon  completing  that  part,  we 
will  turn  our  attention  to  the  main  problem:  identifying  the  largest  integer  k  for  which  the  uniform 
boundedness  of  the  sequences  {mklh}h  is  satisfied. 

Lemma  3.11.  Let  { Mk}h  be  any  box  spline  scale,  and  let  \  be  the  characteristic  function  of  a 
O-neighberhood  B.  Then,  for  sufficiently  small  B,  the  operators 

Jh-  f  ^  (X*i  fhUMhY 

satisfy  the  four  requirements  of  Theorem  2.6. 

Proof.  Jh  certainly  has  the  form  required  in  (2.5),  with  T \  =  x/M/,.  Condition  (b)  there 
can  be  satisfied  by  ensuring  supp\  C  [~x  ..n]d.  Condition  (d)  (listed  in  Theorem  2.6)  holds,  since 
on  B  suppy,  1  -  A//,7\  =  0.  It  remains  to  deal  with  condition  (c)  of  (2.4).  For  that,  we  first 
observe  that  since 

KM  =  TT  f  e(h A< dt, 

{A/*}/,  converges  uniformly  on  suppx  to 

1  dt, 

and  this  latter  expression  is  bounded  away  of  zero  on,  say,  [— ir..Tr)d.  Thus,  for  sufficiently  ho, 
(7),  =  \f Mh}h<h0  are  uniformly  bounded,  as  required.  4 

In  order  to  deal  with  the  essential  requirement  of  Theorem  2.6,  that  is  the  uniform  boundedness 
in  £q(2ir7Zd\0)  of  the  sequences  {mk,h}h,  we  will  prove  the  following: 

Lemma  3.12.  Let  {A/aJa  be  a  box  spline  scale  generated  by  M  —  Mz,\,  a nd  let  Mo  be  the 
corresponding  polynomial  box  spline,  i.e.,  Mo  =  M-=,o-  Given  an  origin-neighborhood  B  C  [— 7r..7r]d, 
define,  as  in  Theorem  2.6, 

mkth  :  (27r2Zd\0)  9  /3  ~  \\(h  +  \-\)kK( •  + 

If,  for  1  <  q  <  oo,  Mq  €  L,(IRrf),  then,  for  small  enough  ho,  the  sequences  { mkth}h<h0  are  uniformly 
bounded  in  f,(2?rZZd\0),  for  k  :=  k'(E). 

Proof.  cf.  §4.2.  4 

With  the  aid  of  the  two  last  lemmas,  we  establish  the  following  theorem: 


9 


Theorem  3.13.  Let  {Mh}h  be  the  box  spline  scale  associated  with  the  box  spline  M  :=  Me, a- 
Consider  the  approximation  maps 

f  «  07 »A*(/), 

with 

Ml)  =  £  W»(-  -  a)J*(/)(o), 

a€2Z<1 

where 

Ml)  ==  (x//A4)\ 

and  \  is  tin*  support  function  of  some  origin-neighborhood  B  C  [— 7r..7r]d.  Then,  for  fc  =  jfc'(E), 
ami  2  <  p  <  oc, 

11/  -  fffc-4h(/)||p  <  constpj|/||piifc/ifc, 

for  every  /  6  V'p,a.  and  every  small  enough  h.  The  result  is  valid  for  p  =  oo,  as  well,  provided  that 
A/o  €  1 1,  with  A/o  =  A/e,o  the  associated  polynomial  box  spline. 

Proof.  The  claim  of  the  present  theorem  with  respect  to  p  =  oo  follows  by  an  application 
of  Lemma  3.12  and  Lemma  3.11  to  Theorem  2.6.  The  same  is  true  also  for  2  <  p  <  oo,  as  soon  as 
we  show  that  for  such  p.  Mo  necessarily  lies  in  Lq.  Since  we  assume  that  the  rank  condition 

rankE  =  d 


holds,  we  can  find  a  d  x  d  invertible  submatrix  X  C  E.  Then 


Afo  =  Mi  *  M2, 


with  Afi  :=  A/.y,o  and  A/2  =  M(£\X), o*  The  trivial  bound 


|  fl  e-*"*  dt\  <  1,  £,u/  G  IRd, 
_  Jo 


proves  that  (*>  G  L 00  fmany  polynomial  box  spline  (f>,  and  therefore  A/2  G  Loo.  Consequently,  it 
suffices  to  prove  that  M\  G  Lq.  By  applying  a  linear  change  of  variables,  we  may  assume  without 


loss  that  X  is  the  identity  matrix,  and  thus  M\  becomes  the  tensor  product  of  the  univariate 
function 


f 1 

uJ  *-*■ 

Jo 


e 


~iwt  dt  = 


1  —  e‘ 


IU1 


which  lies  in  Lq(IR)  for  q  >  1.  It  follows  that  Mi,  hence  also  Mo  are  in  Lg(IRrf). 


* 


As  mentioned  in  the  introduction,  results  similar  to  Theorem  3.4  were  derived  in  [BR2],  but 
with  respect  to  oo-norm.  It  is  shown  there  that  the  Loo-approximation  order  provided  by  M  is 
always  bounded  above  by  fc'(E),  and  this  bound  is  proved  there  to  be  the  exact  approximation 
order  under  the  additional  assumption 

(3.i4)  y.  n 

/362Z<‘\0€€H,{/3?t0 

In  comparison,  Theorem  3.13  requires  the  more  verifiable  condition  Mo  G  L\  (which  is  shown  to 
imply  (3.14)).  For  example,  the  latter  condition  is  satisfied  whenever  E  can  be  partitioned  into  two 
matrices  E  =  X  U  Y  both  of  rank  d.  Indeed,  the  full  rank  assumption  on  X  and  yjmplies  that 
Afx.o,  My-,0  G  L2,  and  hence  their  convolution  product  Ms,o  is  in  the  algebra  A  :=  {/  :  /  G  Li}. 
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Example  3.15.  Let  d  =  2.  Then  one  easily  checks  that  E  can  always  be  partitioned  into  two 
matrices  of  rank  2  unless  E  =  .V  U  {£}  with  .V  is  a  rank-1  matrix,  or,  in  terms  of  Jt(Ei  (cf.  (3.U)i  if 
and  only  if  k lE)  =  1.  As  observed  in  the  proof  of  Corollary  3.8.  fc'(E)  <  k f  E).  hence  in  such  a  case 
k'\E  €  {0. 1}.  Sow.  if  k'(E)  =  0.  then,  since  the  Z  *, -approximation  order  is  proved  in  [BR2J  to 
be  bounded  abate  by  Ar'tE).  we  conclude  ,  tat  M  provides  approximation  order  0.  Therefore,  the 
only  bivariate  box  splines  whose  .^-approximation  order  cannot  be  decided  directly  bv  the  results 
here  are  those  associated  with  a  matrix  E  that  satisfies  k(E)  =  k'(E i  —  l. 


Next,  we  want  to  show,  with  the  aid  of  an  example,  that  the  approximation  order  k'i  can 
exceed  in  times  the  local  approximation  order  of  the  space  H  .U)  fcf.  third  paragraph  of  §3.2!: 

Example  3.16.  Let  d  —  2, 


1/2  0  1/2  1  1/2 

0  1/2  1/2  1/2  l 


1/2  \ 

1/2/ 


and  A  yet  to  be  determined.  It  is  easily  verified  that  Ar'(E)  =  1  <  kiE)  =  o.  and  hence,  by  Theorem 
3.5.  and  in  view  of  the  previous  example,  the  /^-approximation  order  provided  by  is  1.  for 

every  2  <  p  <  x  { and  regardless  of  the  choice  of  A).  The  choice  A  -  0  leads  to  a  stationary 
situation,  and  the  approximation  order  1  must  then  imply  that  the  shifts  of  ,1/r  0  partition  the 
ur.ity,  as  one  can  verify,  with  the  aid  of  Poisson's  summation  formula,  from  the  fact  that 


A/h.  o!2~Ji  =  6. 


J.o- 


J  6  2Z 


On  •lie  other  hand,  for  a  generic  choice  of  A.  H(.!/)  —  (()}.  This  can  be  proved  as  follows,  if  Hi  A/ 
is  turn-trivial,  then  by  Lemma  3.1  of  [BAR],  it.  contains  an  exponential  _  —  e1®  *.  The  frequency 
tO  of  that  exponential  must  satisfy  the  following  two  conditions: 
i a )  A.v  -  i\r9  =  0.  for  some  2  x  2  X  C  E  of  rank  2  ;  and 
;  !>  .Uh.v#  *  2trJ)  =  0.  J  G  ZZ2\().  (cf.  [Rl:§4]  for  (a)  and  [BR1:§2]  for  (b).> 

To  see  that  the  above  la)  and  (b)  can  hold  only  in  exceptional  circumstances,  we  proceed  as  follows: 
we  fix  tin-  above  E.  X  and  A.v-  This  determines  a  unique  ft  (see  (a)  above).  We  now  try  to  define 
y  .  (  6  E,  so  that  |b)  above  is  valid.  For  that,  one  verifies  first  (directly)  that,  regardless  of  the 
specific  choice  of  the  2x2  X  C  E.  M\.\x  can  vanish  only  on  a  subset  of  some  proper  sublattice 
C \  of  2~7Z~.  Selecting  3  G  2~S"’\{0}  in  the  complement  of  this  sublattice,  (b)  above  implies  the 
existence  <>f  <  <=  E\.V  such  that  dt  =  0,  or,  equivalently. 

A(  —  ■  (0  +  3)  G  2?r2Z. 

This  shows  that  one  of  A^,  C  G  E  must,  be  chosen  from  the  countable  set  —i(  ■  (9  +  2t2Z‘),  and 
t.hererfore,  enerically,  H (.If)  =  {()},  as  claimed.  4 

Finally,  we  state  our  result,  concerning  simultaneous  approximation. 

Theorem  3.17.  In  the  notations  of  Theorem  3.  III.  and  for  2  <  p  <  oo, 

11/  -  crtiAh(f) ||iv-  <  constp  r  /t*~rj|/||J).fc, 

for  every  f  G  lp  *,  every  small  enough  h,  and  every  integer  r  <  k.  The  same  holds  for  p  =  cc, 
provided  that  Mu  and  all  its  derivatives  up  to  order  r  inclusive  lie  in  the  Wiener  algebra  A  (or. 
equivalently.  |  •  |rA/0  G  L\.) 

Proof.  cf  §4.4.  A 
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4.  Proofs 


4.1.  Proofs  of  Theorems  2.6  and  2.9 

The  approximation  maps  {.4*}*,  that  we  employ  are  intimately  related  to  those  used  in  [BR2j. 
In  fact,  the  latter,  albeit  a  special  case  of  the  present  schemes,  seem  to  be  their  most  natural  choice. 
In  contrast,  the  error  analysis  of  [BR2]  cannot  be  adopted  here:  that  analysis  makes  use  of  the 
optimal  approximation  to  the  exponential  functions 

ng  :xH-e:'9x,  delR1' 

of  [R2j.  and  synthesize  those  optimal  approximations  on  the  Fourier  domain  to  yield  optima! 
approximation  to  other  smooth  functions.  However,  when  p  <  x  the  above  exponentials  are  not 
in  LP  anv  more,  hence  the  [BR2]  approach  cannot  go  through.  Instead,  wo  use  here  the  following 
identity 

(4.1:  Y  v(-  -  a)g(a)  =  Y,  v*(e.jg). 

*6'i Z*  J62x'/Ld 

which  is  valid  for  anv  compactly  supported  distribution  v  provided  that  g  is  sufficiently  smooth, 
saw  g  €  C^dR'1)  iff.  Theorem  2.6  of  (RS]f.  The  convergence  of  the  right  hand  side  of  '4.1 ,  is  valid 
in  the  topology  of  tempered  distributions  (and  in  most  circumstances,  in  much  stronger  topologies  a 
provided  that  g  and  all  its  derivatives  grow  no  faster  than  polvnomially  at  x. 

Proof  of  Theorem  2.6.  Fix  2  <  p  <  x .  and  let  /  6  \p,k-  We  try  to  estimate  the  error 
(4.2)  11/  -  <ThUh(/))ilp  =  *rf/pIK ihf  -  .-U(/)llp. 

with  the  right  hand  side  in  (4.2)  obtained  from  the  left  hand  side  by  scaling.  Invoking  (4.1)  with 
respect  to  .4 iJf)  (i.e.  with  g  :=  .//,(/)  and  c  :=  Oh),  we  obtain  that 

\Wi/hf  -  .-U(/)!!„  HK/t>/  -<Ph  *  Jhif)  -  Y  ^  *(e.o4(/))i!p 

36  2*Zd\0 

(4.3j 

<\Wi/hi  iip +  n  Y  ^*(cjA(/))Hp- 

P62>rZ'i\0 

We  estimate  each  of  the  two  terms  in  the  last  line  of  (4.3)  with  the  aid  of  the  Hausdorff- Young 
inequality: 

(4.4)  H/llp  <  constpH/H,,  1/p  +  1/?  =  1, 

valid  for  2  <  p  <  oc.  provided  that  f  €  Lq(fRd). 

We  first  estimate  in  the  proposition  below  the  second  term  in  the  second  line  of  (4.3)  For 
later  use.  we  derive  that  estimate  in  a  slightly  more  general  setup  than  needed  here. 

Proposition  4.5.  Let  w  be  some  sequence  defined  on  2~/Z',\0  and  having  (at  most)  polynomial 
growth.  Then,  in  the  notations  of  Theorem  2.6,  anti  under  the  assumptions  there, 

h'l/p II  Y  «'W)<Ph  *  (e;iJh(f))]\p  <  const 

06  2>r2Z'(\0 


Proof.  There  is  nothing  to  prove  in  ease  ||mfc.fcU'i|f„(2TZ‘1\o)  =  Otherwise,  by  (4.4)  it 
suffices  to  bound 

(4.6)  /i‘</p|j(  ^2  u(3)Oh  *  (ejJh(f))r\\Q  =  hd/p\\  55  w(-J)oh Jh(f)(-  + 

J€'2t/Z‘\0  de2rZ’i\0 


(The  justification  for  the  term-by-term  application  of  the  Fourier  transform  is  given  in  the  sequel :. 
Substituting  h~df(-/h)Th  for  Jh(f)  (cf.  (2.5))  we  obtain  from  (4.6)  the  equivalent  expression 

14-7}  h-*/i  HOfc  51  ui-3)f((-  +  3)/h)  Th[-  +  J)\\q. 

J€2tZ'j\0 


Here,  the  infinite  sum  in  the  above  expression  trivially  converges  (since  supp  Th  C  D  c{~~  . .  -]d), 
and  the  limit  is  supported  in  (2rrZZ'i\0)  +  B.  Also,  since  the  weights  {«•(./)}  are  of  polynomial 
growth,  the  convergence  holds  in  5'.  In  retrospect,  this  justifies  the  term-by-term  application  of 
the  Fourier  transform  in  (4.6),  as  well  as  the  changing  of  the  order  of  summation  and  multiplication 
by  Ot,  is  the  display  afterwards. 

We  fix  a  €  2~2Z't\0.  and  compute  that 


i|oh  51  u.'(-.iV((-  +  3)/A)rh(-  +  J)||<f(o+B; 

Je  2rriZ  ‘\o 

=  !ivia)oh/((-  -  a)/h)Th(-  -  q)H‘[^{o+0) 

<  viar^Wl^  ll(A  +  H)*/(-A)ill,(B)  II (h  +  |-|)-*^(-  +  o)\\l^B) 

=  IK  I  +  U)kn\l',B,H)  ll<*  -  +a)IILisv 


Summing  over  all  a  €  2rrZZr\0  and  using  the  uniform  boundedness  of  {Th }  we  obtain  the  bound 
h-'^u  55  u:(-.3)f«,  +  mm-  +  <  const  /t*’H/!l’*!|/«*.fc»r||J4(2!rZ,V0). 

3e_2-7J.‘\o 


and  the  required  result  follows.  4|k 

In  view  of  (4.2),  (4.3)  and  the  claim  of  the  last,  proposition  for  the  choice  u:  —  1,  the  proof  of 
Theorem  2.6  is  reduced  to  the  study  of  the  first  term  on  second  line  of  (4.3).  Here  we  have,  for 
some  positive  const, 

const!|er1/h/  -  <ph  *  Jh(f) ||p  < 


(4-8)  \Wl/hf  -  4>hJh{f)\\q  <  Ikl/h/(l  -  to)||£,(B)  4-  lkl//>/llt,(IR«\B), 

with  the  first  inequality  by  (4.4)  and  the  second  equality  because  Jhif)  —  Th<0\/hj  and  TX  is 
supported  in  B.  Changing  variables,  the  term  l|<H/ft/llz,,(irt''\B)  can  t>e  bounded  as  follows: 

=  h  d  +  'i'/?||/||BQ(lR'<\(B/h)) 

<h~d/p(l  +  c/h)~k ||(1  4-  |.|)*/|| t,(Il^(B/i.))  =  h~d/phk\\f\\p,k  o{l). 
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with  the  o(l)  expression  uniformly  bounded  in  h  and  /  €  Vp.k-  As  for  the  other  term  in  (4.8). 

ikt/h/U  —  °hTh)\\iq(B)  <  |l(ft  +  |'|)*0'l/h/lk,(B)l|(/1  +  H)"k(l  ~  4>hThY\Lv{8»- 

By  assumption  (d),  \\{h  +  )-i)-fc(l  —  O/iTT,  1  [[ l-k(B)  is  bounded  independently  of  h.  while  a  change  of 
variables  yields 

ill*  -  H)*<WI!i,(B)  -  h-dm  +  ttMh'iB)  =  **-,i/p!Ki  +  H)7;k.,sfi/h,  < 

Thus,  we  have  shown  that  each  of  the  two  terms  in  the  last  line  of  (4.3)  is  of  order  0{hk~d-P).  and 
the  claim  of  Theorem  2.6  then  follows.  A 

Proof  of  Theorem  2.9.  The  proof  of  Theorem  2.9  closely  follows  that  of  its  special  case.  The¬ 
orem  2.6.  We  thus  only  outline  the  proof,  emphasizing  parts  of  the  proof  that  deviate  from  their 
counterparts  iri  Theorem  2.6. 

We  let  P  €  Ilr  be  a  homogeneous  polynomial,  and  P[D)  being  its  associated  constant - 
coefficient  differential  operator.  We  want  to  establish  the  bound 

i.4.0i  \\P{D)(f  -  <7h.4h(/))||p  <  const  hk~dr* 7/WiK.J- 

Such  an  estimate  leads  to  the  desired  result,  since  we  may  range  P  over  some  fixed  homogeneous 
basis  for  Pr.  Here,  const  should  be  independent  of  h  and  /.  but  may  depend  on  P.  p.  and  r. 

For  the  proof  of  (4.9).  we  consider 

P(D)(f  -<xhAh(f))  =  P(D)f  -h-At'pohP(D)AH{f). 

By  (4.1)  (with  c  :=  <f>h  and  g  :=  Jh{J )), 

(4.10)  P(D)Ah(f)  ~  P{D)  Y,  4>k*(ejJh(f))  =  Y  ^*(p( D)(ejJh(f))). 

ZZl-71:1  3£2x?/.i 

with  the  changing  of  order  of  summation  and  differentiation  justified  by  the  ^'-convergence  of  the 
sum.  The  term  corresponding  to  3  =  0  is  <ph  *  ( P{D).Jfl{f ))  =  /i<leg  p0il  *  .Jh(P(D)f),  and  we  first 
estimate 

P{D)f  -  Oh{<t>h  *  Jh(P(D)f)). 

Since  /  £  V’p*,  P{D)f  €  V’p,it-deg  p.  and  hence  the  proof  of  Theorem  2.6  yields  that  for  small 
enough  h  we  have 


II P(D)f  -  oh(4>h  »  MP(D)f)) ||p  $  const/tfc-de*p||P(£2)/||p.fc_ 


p.fc  —  deg  P- 


It.  remains  t.o  bound  the  expression 


Y  <Ph*P{D)(eaJh{f)). 

BZ2*T7P\<i 


Here,  we  fix  fl  ^  0,  and  expand 


P(D)^Jh(f))  =  Y,  ~y  (D"P)iD)(e0)  D"Jh(f)  =  e Y  ^  (0°P)(i.J)  DaJH(f), 


and  use  summation  by  parts  (allowed,  since  the  range  of  a  is  actually  finite,  and  since,  for  any  fixed 
a,  the  Fourier  transform  of  the  various  ,3-summands  have  pairwise  disjoint  supports)  followed  by 
the  triangle  inequality  to  estimate  the  norm  of  that  part  as  follows: 

n  £  a  •  p(DX'<,j*(m, 

0G2w2Z*\O 

(4.H)  *£^11  £ 

a>0  ‘  062*2^0 

=£^n  £  (.D-p)mA>M<.D”fmt 

a>0  '  0£2ic2Zd\O 


Invoking  Proposition  4.5  with  w(0)  :=  wa(0)  :=  ( DaP)(i0 )  (and  with  /,  k  replaced  by  Daf  and 
k  —  |a|i  respectively),  we  obtain  that 


(4.12) 


A'“l,ll  £  (£>“P)(i/3)^»(e«A(C”/))llf  = 

/3€2jr2Zd\0 

*W‘!I  X  wa(P)<t>h  *  (e0Jh(Da f))\\p  <consthfc||Da/||P(fc_|a|J|mfc,hu;a||<<(2,z^0). 

0&*7L*\O 


Since  (i):  the  actual  range  of  a  in  (4.11)  is  finite  (|a|i  <  degP),  (ii):  || JZ>“/||J,,*r_|0|l  <  ||/||Plfc,  and 
(iii):  \(DaP)(i0)\  <  const|/3|r,  we  derive  from  (4.11)  and  (4.12)  the  inequality 


||  £  A*  P(D)MUmp  <  cOQStA‘||/||„,1||m;.J|,)(2.z.w, 

/?G2ir2Z‘,\0 

from  which  (4.9)  follows. 
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4.2.  Proof  of  Lemma  3.12 

Given  0  G  2-rr7Zd\0,  our  first  goal  is  to  estimate  |M/l(u/  +  /?)|,  for  small  u.  Initially,  this  is  done 
without  using  the  assumption  Mq  €  Lq.  We  consider,  one  by  one,  the  factors 

(4.13)  |  f  dt\ 

Jo 

that  form  |Af*|.  For  this,  fixing  0  G  27r2Zrf\0,  we  partition  5  according  to  the  behaviour  of  their 
corresponding  factors  into  three  groups:  the  first  of  which  is  Kp  (cf.  (3.6))  and  the  other  two  are 
defined  as  follows 

£„:={£€  S:  £-0*2Z}> 

Op~{(€Z:  £•/?  =  ()}. 

Case  I :  £  €  Op.  In  this  case,  for  sufficiently  small  h,  we  have 

|  f  c(A>e-*<  (w+^))‘  dt\  <  2. 

Jo 
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Case  II:  f  €  Lj.  Here  we  use  the  estimate 


(4.14) 


J  f  e(hA < 
JO 


_ 3 _ 

|AA?  -  if  •  (jj  +  J)| 


valid  for  all  *  and  sufficiently  small  h.  The  rationality  of  E  implies  the  existence  of  a  €  TL  such 
that  nE  is  integral,  and  thus  we  have  If  •  J|  >  2 ~/n,  and  (4.14)  shows  that  for  sufficiently  small  h 
and  v 


If  ■  J| ' 


with  c  depending  only  on  r;. 

Case  III:  f  €  Kj.  In  this  final  case,  we  write 


(4.15) 


t  ' 


dt\  = 


_  jj 

I h\  -  if  •  (w  +•  J)| 


|ehAe-1^*'  -  11 
\h\t  -  if  •  -  J)j 


The  denominator  in  the  right  hand  side  of  (4.15)  can  be  estimated  as  in  the  previous  case,  while 
the  numerator,  for  sufficiently  small  h  and  jw'|.  can  be  bounded  by  c(h  +  ;^j).  hence  we  obtain  in 
this  case  the  estimate 


j  df{  <  c> 

Jo 


(h  +  M) 

If  '  Jl 


Combining  these  various  estimates  we  obtain  that,  for  uj  in  some  h-independent  neighborhood 
of  the  origin,  for  h  sufficiently  small  and  for  some  (/t. j)-independent  const. 


(4.16)  |A/^(u;  4-  (3)|  <  const(/t  +  |f  •  0\~ 1  <  constf/i  +  |w!)fc'(S>  ][[  !f-i!_i. 

£  €  K j u  Lj  £€  K ju L  j 


I 


I 


I 


the  second  inequality  since  k'(-)  <  #I\j  for  every  3  €  2~ZHd\0. 

For  later  use,  we  record  this  intermediate  estimate: 

Proposition  4.17.  Let  (A/),  f/j  fee  a  box  spline  scale  associated  with  E  and  A.  Then ,  for  k  :=  k'{E). 
for  some  origin-neighborhood  B,  for  sufficiently  small  h,  and  for  some  ( h .  3)-independen t  const  we 
have 

mk.h(3)  :=  ||(/i  +  |.|r*A/a(.+/3)||tM(B)  <  const  ]J  |f  -  ^|~l. 

£  €  K  j  u  L  j 

To  complete  the  proof  of  Lemma  3.12,  we  need  to  show  that  the  sequence 
<h?  =  n  £e27r2Zd\0 

K jU Lq 

is  in  f,(2~/Zd\0),  whenever  M0  €  Lq.  A  slightly  more  general  assertion  is  proved  in  the  following 
proposition. 

Proposition  4.18.  Assume  that,  for  the  polynomial  box  spline  A/o  and  for  some  r  >  0.  |-|r-U»  € 
Lq.  Then 


(4.19) 


£  mrq  n  if-^r?<oc. 

0eM/.d\O  (€Kj,u!.a 


» 
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Proof  of  the  Proposition.  Choose  u  €  IRrf  that  satisfies  the  following  conditions: 

(a)  For  every  £  6  E,  4  •  w  is  not  27T-rational. 

(b)  The  series 

(4.20)  Y1  iriwo(w+u)r 

converges. 

(c) 

K-(u/  +  /?)|<2|£-£|,  £  G  E,  0€2ir7Zd, 

It  is  clear  that  the  set  of  points  u>  €  IRrf  that  violate  (a)  is  of  measure  0.  Also,  because 
| •  | r A/o  €  L,(IR'<),  the  series 

Y,  iK+^nvo^+w 

0€2nZd\O 

converges  a.e.  on  [— n..7r]d~  implying  thereby  the  a.e.  convergence  of  (4.20).  On  the  other  hand, 
because  £  6  E  is  rational,  inf{|£  ■  0\  :  0  €  27r2Zrf\0,  £  •  0  ^  0}  >  0,  and  hence  condition  (c)  is 
satisfied  by  all  small  enough  u.  This  proves  that  there  is  u  that  satisfies  all  the  above  conditions. 
The  rationality  of  E  implies  that  for  any  fixed  £  €  E  the  range  of  the  map 

(27rZZrf\0)  9  0  »-*  e~iH“+0)  -  1 


is  finite,  and,  because  of  condition  (a),  this  range  does  not  contain  0.  Condition  (a)  also  implies 
that  fQl  e~i(  u,tdt  ±  0,  for  all  £  €  E.  Thus,  we  obtain  from  condition  (c)  the  estimate 


(4.21) 


n 

teKgULg 


g -i£  (ul+0)  _  i 

f  •  (a/  +  0) 


>  const  |£  •  0\ 

teKgULg 


Condition  (b)  then  implies  the  desired  result. 
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4.3.  Proofs  of  Theorems  3.4  and  3.7 

The  positive  statement  in  Theorem  3.4,  i.e.,  that  the  approximation  order  provided  by  Mh,a 
(to  functions  in  V2tk^s)  —  W”*  is  at  * east  *'(—)»  follows  from  Theorem  2.6  when  combined  with 

Lemma  3.11  and  Lemma  3.12.  Indeed,  we  only  need  to  verify  that  the  requirement  M0  €  L2  (needed 
for  the  application  of  Lemma  3.12)  holds.  That  was  proved  in  Theorem  3.13,  but,  as  matter  of 
fact,  also  follows  directly  from  the  fact  that  Mo  €  L 2. 

The  negative  statement  of  Theorem  3.4  will  follow  from  Theorem  3.7.  Indeed,  Theorem  3.7 
provides  jfc'(E)  as  an  upper  bound  on  the  approximation  order  (and  in  the  strongest  possible  sense). 

Therefore,  only  Theorem  3.7  requires  a  proof. 

For  that  proof,  we  need  to  borrow  some  of  the  general  tools  and  results  developed  in  [BDR1]. 
We  remark  that  the  definition  of  the  PSI  space  S(<f>)  in  [BDR1]  differs  from  the  one  given  here: 
it  is  defined  there  as  the  /^-closure  of  the  algebraic  span  of  the  shifts  of  (f>.  However,  as  Theorem 
2.13  of  [BDR1]  asserts,  S(<f>)  T\L2  of  the  present  paper  is  dense  in  S(<j>)  of  [BDR1],  hence  the  two 
spaces  share  the  same  approximation  orders,  and  to  the  same  functions. 
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We  first  define  for  every  h  >  0 


.\k  :=  A4*A/h(~). 

Note  that 

Xh  =  |A4|2- 

Because  of  the  rank  assumption  (3.2).  A£,  6  £2,  and  therefore  AJ  (which  is  clearly  compactly 
supported)  is  continuous,  as  any  convolution  product  of  two  £  >-functions  is.  We  make  a  substantial 
use  of  the  symbol  A7,  of  A'/,  defined  as 

Xh  :=  Y.  'Vh(a)e_0. 

a£TLl 

A  standard  application  of  Poisson's  summation  formula  shows  that 

-v*=  Y  - •*)  =  E  IA4('  +  ^)!2- 

3€2z7Z'‘  3€2*7ld 

Thus,  we  see  that  Xh  is  a  non-negative  trigonometric  polynomial,  and.  further.  .Y h{y)  =  Q  only  if 
Xk(ij)  =  0. 

The  £.>-approxiination  orders  provided  by  M  are  determined  [BDR1],  by  the  behaviour  around 
the  origin  of  the  functions 

(4.22)  Ah:=(l-Xh/.\)l/\  h>  0 

(here,  0/0  is  defined  as  zero,  but  in  any  case.  X\,  as  a  trigonometric  polynomial,  vanishes  on 
a  null-set).  Note  that  \k  is  non-negative  and  bounded  by  1.  The  precise  result  that  we  need  here 
follows  from  Theorem  2.20  and  Corollary  3.10  of  [BDP1],  and  reads  as  follows: 

Result  4.23.  Let  f  €  W*  •  Let  {<?h}h  be  a  subset  of  Li-  Then 
(4.24)  dist2  (f^k(S(d>h)))=o(hk) 

onlv  if 

h-d/2ll‘U*Jh,(B)  =  o(hk), 

on  some  or igin- neighborhood  B. 

We  next  attempt  Jo  replace  A/,  in  this  result  by  simpler  expressions.  It  is  clear  that  we  might 
replace  A/,  by  (AJ  -  AJ)1/2  in  case  the  ratio 

(Nh-Nh)l/2  zrU2 

- - - —  =  At, 

Ah 

is  bounded  around  the  origin  by  h-independent  positive  constants.  For  this  we  need  the  following 
lemma,  in  which  we  make  use  of  the  fact  that  every  A J,  and  in  particular  .Vp.  is  supported  in  the 
symmetric  region  Z^^_^  —  Z-z  -  Z=  (which  follows  from  the  fact  that  A£(— -)  =  cf.  (3.3) 

for  the  definition  of  Z-=). 


Lemma  4.25.  Nh  -+  N '0  uniformly  and  hence  Nh  -*•  iVo  in  any  p-norm,  1  <  p  <  oo. 

Proof.  The  first  claim  easily  follows  from  the  distributional  definition  of  box  splines  (cf. 
e.g.  Definition  2.1  in  [Rl],  More  specifically,  one  can  apply  to  Nh  -  N0  the  argument  used  in  the 
proof  of  Lemma  5.1  in  [DR]).  The  second  claim  follows  from  the  first,  since  all  Nh  are  supported 
in  the  same  compact  domain  Z-=  -  Z-=.  4 

Thanks  to  Lemma  4.25,  we  know  that  {Nh.}h  are  uniformly  bounded  around  the  origin,  hence 
may  replace  {A*}/,  in  Result  4.23  by 

v*)"2  =  (  Y,  |AW-  +  /3)I2)1/2. 

^627rZ'1\0 

Thus,  if,  for  some  /  e  Wk ,  dist2(/,  <Xfc(S(</>/i)))  =  o(hk) ,  then  we  must  have,  for  every  0  6  27rZZrf\0, 

+  =«<**). 


which  implies  by  scaling  that 

(4.26)  ||AWft  ■  +/3)/IUl(A-.e,  =  o(h“). 


Let  0  e  2a-ZZrf\0  be  chosen  with  #K0  =  k.  Since  /  6  L^,  it  is  supported  on  a  set  of  positive 
measure,  and  therefore,  for  sufficiently  small  £,  the  set 

Ae  :=  {uj  €  IR*  :  |A€  -  % -u\  >  e,  €  5} 

has  a  positive  measure  intersection  with  supp /.  We  fix  such  e,  and  we  let  Q  be  any  bounded 
measurable  subset  of  Ae  for  which  ||/||z.2(rj)  >  0. 

A  straightforward  computation  shows  that,  for  any  9  €  IR5  and  £  €  E\Ap(E), 


(4.27)  [' dt  — „  f 

Jo  /l-,°  Jo 

uniformly.  Further,  for  £  €  K0,  we  get  that 


v: 


e(h\(-n-(hff+0))t  dt 


1  A;  -  •  9 


>  const  >0,  V0  6  Ac, 


a-*o  I  £  •  0 

wkh  the  convergence  being  uniform  on  compact  sets,  hence  on  Q.  Thus,  for  small  enough  ho, 
m{{h~l\£  e(hx(-iHh0+0»tdt\  :  h<  h0}  >0,  €  K0. 


This,  together  with  (4.27),  implies,  with 


rh  :=  inf{|M/,(h0  +  /?)]:  9  €  Q}, 


that,  for  small  enough  h, 


h~krh 


9  6  ft}  >  const  >  0. 


19 


Thus,  by  (4.26). 

o(kk)  =  p/*(A  ■  +J)/||i,(n)  >  ll/llij(n)rA  >  const||/||£,(n)/i*. 
a  contradiction  to  the  fact  that  ||/!li,(ni  /  0. 

This  completes  the  proof  of  Theorem  3.7.  and  thereby  the  proof  of  Theorem  3.4.  4 

4.4.  Proof  of  Theorem  3.17 

The  proof  invokes  Theorem  2.9  for  the  choice  Oh  Obi ■  Thus,  we  need  to  verify  that  for  an 
integer  r  <  k  (with  k  Ar'(H)).  the  sequences 

m[  h  :  (2-ZZd\0)  3  3  -  \3\r\\(h  +  h  <  hQ 

lie  in  /,j( 2t2Z'i\0)  and  are  bounded  there. 

For  that,  we  first  invoke  Proposition  4.17  to  conclude  that 

mrk  ^(3)  <  constl J|r  JJ  •  d|_1. 

A'juL 

Therefore,  by  Proposition  4.18.  the  uniform  boundedness  of  {mrkh}k  in  (,li'2~7Z,l\Q)  is  implied  by 

the  condition  Mr.\/o  €  L,r  For  p  -  oc  (i.e..  q  =  l)  this  latter  condition  is  assumed  in  the  present 
theorem.  \\>  prove  here  the  validi'v  of  the  condition  for  n  >  1. 

We  first  observe  that  the  condition  the  condition  [•|r.l/o  is  equivalent  to  the  statement:  "for 
each  homogeneous  polynomial  P  of  degree  r.  (P{D)MoP £  Lr"  Here,  r  <  fe'iH)  -  1  <  k (E)  -  1 
icf.  the  proof  of  Corollary  3.8).  On  the  other  hand,  it  is  known.  [BH],  that  for  any  polynomial  P 
of  degree  <  k( E),  P(D)Mq  can  be  written  as  a  finite  sum 

(4.28)  ^c.v.Ux(-  -  q.y), 

x 

where  each  A'  in  the  above  sum  is  a  submatrix  of  E  of  full  rank  d.  where  c\  are  some  coefficients, 
o.v  €  1R  ,  and  M\  =  Mx.o  is  the  polynomial  box  spline  defined  by  A*.  Since  each  X  above 
satisfies  rank  A"  =  d,  then,  as  established  in  the  proof  of  Theorem  3.13,  M\  €  Lq  for  every 
q  >  1.  Therefore,  the  Fourier  transform  of  the  sum  in  (4.28)  is  in  Lq ,  namely.  {P(D)Mq]~€  Lq. 
Consequently,  |- W.V/o  €  Lq.  A 
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